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Within a multiphase transport model we study phi meson production in relativistic heavy ion 
collisions from both superposition of initial multiple proton-proton interactions and the secondary 
collisions in the produced hadronic matter. The yield of phi mesons is then reconstructed from their 
decaying product of either the kaon-antikaon pairs or the dimuon pairs. Since the kaon-antikaon 
pairs at midrapidity with low transverse momenta are predominantly rescattered or absorbed in the 
hadronic medium, they can not be used to reconstruct the phi meson and lead thus to a smaller 
reconstructed phi meson yield than that reconstructed from the dimuon channel. With in-medium 
mass modifications of kaons and (f) mesons, the (f) yield from dimuons is further enhanced compared 
to that from the kaon-antikaon pairs. The model result is compared with the experimental data 
f\J , at the CERN/SPS and RHIC energies and its implications to quark-gluon plasma formation are 

■ discussed. 
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I. INTRODUCTION 



One of the major goals in relativistic heavy ion collisions is to study the properties of hot and dense matter 
and possibly to create and identify a new form of matter, the quark gluon plasma (QGP). It was suggested |^ 
' that enhanced strangeness production could serve as an important signal for the deconfined matter. The dominant 
production of ss pairs via gluon-gluon interaction may lead to strangeness (chemical and flavor) equilibration times 
comparable to the lifetime of the plasma and much shorter than that of a thermally equilibrated hadronic fireball. The 
, subsequent hadronization is then expected to result in an enhanced production of strange and multistrange particles 
and antiparticles. In particular, it has been argued that with the formation of quark-gluon plasma not only the 
production of phi mesons, which consist of ss, is enhanced but they also retain the information on the conditions of 
the hot plasma as it is believed that phi mesons interact weakly in the hadronic matter and therefore freeze out quite 
[ early from the system Q. 

^ , Using a proton or sulfur projectile at 200 GeV/nucleon and a tungsten or uranium target, phi meson production 
^ ■ has been previously studied at CERN/SPS by the NA38 Collaboration and the IIELIOS-3 Collaboration via 
H the dimuon invariant mass spectra. The double ratio + p^))su(W) /{4>l + P^))pW has been measured and was 

^ found to have a value of 2 ^ 3. Various theoretical attempts have been made to understand this enhancement. In 
\ particular, an enhancement of phi meson yield may be a signature of the formation of a quark gluon plasma in the 
1^' > collisions [^. However, this enhancement can also be explained either using hadronic models if one takes into account 
the reduced phi meson mass in nuclear medium |^ or via the fragmentation of color ropes that are formed from the 
initial strings 

^ ' Recently, phi meson production has also been measured in central Pb-|-Pb collisions at 158A GeV at the CERN/SPS. 
5^ , The NA49 Collaboration Q has identified the phi meson via the decay channel </> K~^K~, while the NA50 
Collaboration ||^ measured it using (j) — > fJ-^ fJ-^ decay. It was found that the extracted number of phi meson from 
dimuon channel exceeds by a factor between two and four from that extracted from the K~ channel. This difference 
has been attributed to the fact that not all phi mesons can be reconstructed from K~^K~ resulting from their decays 
due to rescattering of the and K~ in the hadronic matter. The suppression factor can be 40-60% based on the 
RQMD model ^ . This effect has been further studied in a schematic model by including the effect of changing kaon 
masses in the hadronic medium [Q. 

In the present paper, using a multiphase transport (AMPT) model we make a detailed and consistent study of the 
production of phi mesons reconstructed from the K~^K^ and fJ-^ fJ-~ decay channels. The effects of in-medium mass 
modification to the phi meson yield and spectra are also studied and are compared with the experimental findings 
by the NA49 and NA50 Collaborations at the SPS energies. We also present the phi meson yield and spectra at the 
RHIC energy of = 130 A GeV. 

This paper is organized as follows. In Section |l|, we briefly review the AMPT model, discuss the phi meson 
production and interaction cross sections in hadronic matter as well as the medium effects on kaon and phi meson 
masses. The identification of phi mesons from both the dimuon and dikaon channels in the transport model is described 



in Section HI. In section IV, we show the results on phi mesons from heavy ion collisions at both SPS and RHIC. 
Finally, summary and conclusions are given in Section 0. 
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II. THE MODEL 



A. A Multiphase transport model 

To describe the heavy ion coUision dynamics, we employ a multiphase transport model (AMPT) that includes both 
the initial partonic and final hadronic interactions. The AMPT model is a hybrid model that uses as input 
the minijet partons from the hard processes and the strings from the soft processes in the HIJING model The 
dynamical evolution of partons are then modeled by the ZPC jl^ parton cascade model, while the transition from 
the partonic matter to the hadronic matter is based on the Lund string fragmentation model p^ . The final-state 
hadronic scatterings are modeled by the ART model ||l5|. The AMPT model has been quite successful in describing 
the measured transverse momenta of pions and kaons, and the rapidity distributions of charge particles as well 

as the particle to antiparticle ratios |1^ in heavy ion collisions at both SPS and RHIC. Including melting of initial 
strings to partons, the model can also account for the observed large elliptic flow figj ] and the measured two-pion 
correlation function at RHIC. In the original AMPT model, phi mesons are produced only from initial string 
fragmentation. In the present work, we extend it to also include phi meson production and interactions in the hadronic 
matter. 



B. phi meson production and interactions 



For phi meson production from hadron scatterings in the ART model, we consider both the baryon-baryon in- 
teraction channels BB (j)NN and the meson-baryon channels (tt, p)B (j>B, where B = N,A,N*. The cross 



sections for these processes have been evaluated in the one-boson-exchange model |21| and are used here. We have 
also included phi meson production from the reaction KA (f>N with the cross section taken from Ref. |^ based on 
a kaon-exchange model. 

Phi meson can be further produced from kaon-antikaon scattering. For a kaon-antikaon pair with invariant mass 
M, the total cross section for KK — > i/i is taken to be a Breit-Wigner form 

. . ,^^ 37r {m^T^f 

where the phi meson decay width to KK is given by 

/ \ 3/2 

with the coupling constant 9'^j^j^/4:Tt « 1.69 determined from the empirical width of 3.7 MeV, corresponding to 83% 
of the total width, at M = m^. 

For phi meson interactions with baryons, we include the absorption reactions given by the inverse reactions of phi 
meson production from meson-baryon interactions given above. The corresponding cross sections are obtained using 
the detailed balance relations. The cross section for phi meson elastic scattering with a nucleon is taken to be 0.56 
mb as extracted from the data on phi- meson photoproduction using the vector meson dominance model p^ . These 
cross sections are shown in the upper left panel of Fig. 0. 

Phi mesons can also be scattered by mesons. Using effective hadronic Lagrangians, where coupling constants were 
determined from experimental partial decay rates, the total collisional width of cj) due to the reactions 07r KK* , 
(j)p KK, (j)K 4>K and (jxj) KK was found to be less than 35 MeV |Q (where K denotes either a kaon or an 
antikaon as appropriate). However, recent calculations based on a Hidden Local Symmetry Lagrangian shows that 
the collisional rates of (j) with pseudoscalar (tt, K) and vector (p, w, K* , (p) mesons are appreciably large, especially 
for the K* , resulting in a much smaller mean free path, of about 2.4 fm in a hadronic matter at temperature T > 170 
MeV, compared to the typical hadronic system size of ~ 10 — 15 fm created in heavy ion collisions. We have included 
ah these possible interactions, i.e., (/)M -> {K, K*){K, K*), and (f>{K,K*) M{K,K*), where M = {■K,p,uj) with 
cross sections determined from the partial collisional widths given in Ref. HJ] . Specifically, we take matrix elements 
I Alp of these reactions to be independent of the center-of-mass energy and extract their values from the collisional 
widths. The cross sections for these reactions are then given by da/dfl ^ |A^pfc//(fciS^), where s is the cm. energy 
and ki and kf are, respectively, the initial and final cm. momenta of the colliding hadrons. The cross sections for 
the inverse reactions of the above processes are then obtained from the detailed balance relations. In Fig. nl we 
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show the cross sections for phi-pion and phi-rho scattering (upper right panel), phi-kaon scattering (lower left panel), 
and phi-K* scattering (lower right panel). The elastic cross section of phi meson with other mesons can be similarly 
obtained and is found to be about 2 mb. 
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FIG. 1. Phi meson scattering cross sections by nucleons (upper left panel), mesons (upper right panel), kaons (lower left 
panel), and K* (lower right panel). 



C. medium effects 



Since the phi meson mass is only 32 MeV above the kaon-antikaon threshold, it is strongly coupled to the kaon 
and antikaon dynamics, i.e., to the strangeness content of the surrounding medium. A small change in the masses of 
either the phi meson or kaon thus would appreciably affect the decay width of phi meson. 

Studies based on various approaches, such as the relativistic mean- field theory |2^ and the chiral effective Lagrangian 
p6[ , have led to a general consensus that K'^ feels a rather weak repulsive potential while the K~ feels a relatively 
stronger attractive potential. With such modifications, the kaon yield and spectra could be successfully explained at 



SIS energies (1-2 A GeV) p% 



29(1 and even at AGS energy II A GeV) pO]. In the present study, we adopt the kaon 



and antikaon energies from chiral effective Lagrangian, i.e., 

1/2 

where bx = 3/^/8 « 0.333 GeV fm^, ax ~ 0.22 GeV^ fm^ for kaons, and w 0.45 GeV^ fm^ for antikaons. The 
potential, defined by Uk = (^k — \/ "^/f + and given by the plus sign in Eq. (H), is then ~ -1-20 MeV at 
normal nuclear density and is consistent with that determined from the impulsive approximation using the empirical 
kaon-nucleon scattering length The potential, given by the minus sign in Eq. (^), at normal nuclear density is 
about —110 MeV and is somewhat less than that extracted from the kaonic atoms The present set of parameters 
was found to reproduce well the and K~ kinetic energy spectra and the ratio of their yields for heavy ion collisions 
at the SIS energy ||. 
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For the in- medium phi meson mass, the QCD sum rule studies show that it decreases shghtly in hot dense matter 
P^-^ . Using the result from Ref. |Q , we have 

^ w 1 - 0.0255-^ , (4) 

where po is the normal nuclear matter density. 

The in-medium decay width of phi meson to kaon-antikaon pair is then given by 

n^KRiM) = [{{mf - {m*^ + m*^f) {(rn,f - (mj, - . (5) 

Although the phi meson mass decreases in the medium, the larger reduction of the overall kaon-antikaon mass with 
increasing density results in an increase of phi meson width T*^^^ in the medium. At density p — 2pQ, the phi 
meson in-medium width is about 45 MeV. 



III. RECONSTRUCTION OF PHI MESONS 



Since phi meson is unstable, it can only be detected from its decay product of either the kaon-antikaon pair or the 
lepton pair. The kaon-antikaon pair decaying from a phi meson is, however, likely to undergo appreciable rescattering 
in the medium, and this would lead to a reconstructed invariant mass situated outside the original phi meson peak. 
In fact, the momentum transfer to a kaon suffering one collision with a pion at T = 150 MeV can be estimated to be 
about 45 MeV which is much larger than the total phi meson width of = 4.45 MeV. Hence phi mesons decaying 
in the dense medium is difficult to be identified via reconstructed kaon-antikaon pairs. In contrast, dileptons have 
negligible final-state interactions with the surrounding hadronic medium and therefore escape essentially unscathed 
during the entire evolution of the system. D ilep tons are thus considered to carry useful information about hadron 
properties in hot and dense hadronic matter ||37|], which are expected to be different from those in free space. 



A. dileptons 

To detect a phi meson from its dimuon decay, we need the decay width for M^/^^ p^], which in the medium is 

The coefficient C^+p- = a2/27(.g2^^/47r) = 1.634 x 10"^ is determined fr om the measured width at M = m^. Since 
dimuons are emitted continuously during the evolution of the system, the total number of dimuon is given by p8[ | 

^M+M- = iv^(t)r;^^+^- (M) + N^tf) ^^-1'-^^- , (7) 

Jo V p p 

where Nfj^it) denotes the number of phi meson at time t. In the above, the first term corresponds to dimuon production 
before the freeze-out time of tf — 35 fm/c considered in our calculation, while the second term refers to dimuon 
emission after freeze-out. The reconstructed phi meson number is obtained by dividing the above expression by the 
dimuon branching ratio in free-space of F^^^+^-ZF^ — 3.7 x lO^''. Since most previous studies j|Jl^ have neglected 
the phi meson annihilation and production channels via baryons and especially by pseudoscalar and vector mesons, 
the phi meson lifetime in these studies is thus comparable to the lifetime of the system. Consequently, phi mesons 
from dimuon channel get significant contributions only after freeze-out. 



B. kaons 



The number of kaon-antikaon pair stemming from phi meson decay can be similarly expressed as Eq. for dimuon 
production, i.e.. 
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Nkk - r dt N^{t)T* j^{M) + N^{tf)^^t^ . (8) 

The phi meson abundance from kaon-antikaon decay is obtained by dividing Eq. by the KK branching ratio 
in free-space F^^^^/f/F^. It is evident from Eqs. (j^) and (^) that in absence of any medium effect on the masses 
and decay widths, i.e. for F*_^^^ — ^ (f,^KK ^^'^ ^*ci>-~> n+ ^i- ~ 1^^^^+^- , the number of reconstructed phi meson 
from dimuon channel is same as that from the kaon-antikaon channel if all the kaon pairs escape the collision zone 
unscattered. On the other hand, with large in-medium phi meson width both the production and annihilation of phi 
meson from KK is enhanced. If the production of phi meson dominates over its decay especially at the early 
dense stage of the collision, the dimuons originating from the phi meson decay will be enhanced. On the other hand, 
the kaon-antikaon pairs from phi meson decays will be essentially undetected in the phi meson reconstruction as they 
are expected to be rescattered and experience strong medium modifications of their masses. 



IV. RESULTS 



We have used the extended AMPT model described in the above to study phi meson production in heavy ion 
collisions at both SPS and RHIC energies. 



A. SPS 



1. time evolution 



The time evolution of phi meson abundance at midrapidity for central Pb-(-Pb collisions at 158^ GeV in the AMPT 
model is shown in the upper left panel of Fig. |[ The solid curves are the results from calculations where both 
kaon and phi meson medium effects are neglected. The minijet partons and predominantly string decays in the 
HIJING model at the SPS energy produces about 2.6 initial number of phi mesons per event. Since these phi mesons 
are produced at different times as a result of finite formation time, they appear successively in the ART hadronic 
transport. Subsequent multiple rescattering in the hadronic medium leads to a peak value of « 4.0 at f = 3 fm/c, 
which finally decays to ^ 2.3 at an average freeze-out time of i/ ~ 35 fm/c. Because of small decay width of (/) into 
KK pair, most of the phi mesons decay outside the fireball (see upper right panel of Fig. In spite of this small 
width, the production of phi meson from the inverse reaction dominates over the decay due to considerable abundance 
of kaons from the HIJING model at the early stage of the reaction. At times t > 15 fm/c the process </) ^ KK is seen 
to approach chemical equilibrium. Note that even at late times at around the freeze-out value there is a substantial 
phi meson production from the kaon-antikaon channel which has been neglected in the previous study |10|. 

Time evolution of the production and absorption rates of phi meson from different channels are shown in the 
lower panels of Fig. ^j. Of all the coUisional channels considered for the phi meson, the dominant ones are (j)M 
{K, K*){K, K*) and (j)K* M{K,K*) {M = 7r,p, cj), which are of comparable magnitude. Although the cross 
section in the former reaction is smaller compared to that oi (p — K* (see Fig. 1), the larger abundance of pions 
relative to kaons {K'^/ir'^ = 0.12 from the HIJING) leads to similar contribution to phi meson production in both 
these reactions. Moreover, the small abundance of phi meson results in identical production and annihilation rates 
for (f) — M collision (lower left panel), i.e. chemical equilibrium in this reaction is rapidly achieved. In contrast, phi 
meson production overwhelms its destruction (lower right panel) due to larger abundances of the mesons M{= tt, p, uj) 
relative to K and especially K*. In particular, the process MK* (pK* with its largest cross section and smaller 
threshold contributes ~ 80% to the total (p + K* formation. This enhanced production at about 3 fm/c is reflected 
in the rapid increase of the phi meson abundance (upper left panel). The production and annihilation rates for phi 
mesons in processes involving K and K*, which are not shown in the figure, exhibit similar qualitative features as 
for those with K and K*, but are about a factor of 1.5 smaller reflecting the K/K ratio of 0.70 at the SPS energy. 
The phi meson collisional rate as observed here is considerably larger than that predicted in all previous calculations 
1^,^. As the system expands the collisional rate drops and these reactions reach chemical equilibration at t = 7 
fm/c. Finally, at a later time phi meson decay becomes more effective than its absorption. It may be noted that 
the initial dominance of phi meson production and its subsequent absorption as observed in the present transport 
calculation is different from the findings of Ref. [Q, where mesons embedded in a gas of (tt, K, p, uj, K*) at 
chemical equilibrium, are found to decrease with evolution of the system. 

The chemical equilibration of phi mesons via production and annihilation in the hadronic matter can be better seen 
by switching off the initial phi production from string decays in the HIJING model. As shown by the dash-dotted 
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FIG. 2. Upper left panel: Time evolution of midrapidity (|j/| < 0.5) phi meson from Pb+Pb collisions at the SPS energy 
of 15SA GeV at an impact parameter of fe < 3.5 fm in the AMPT model. The results are for without (solid curve) and 
with (dashed curve) in-medium mass modifications, and with further increase of phi meson number by two in the HIJING 
model (dotted curve). The phi meson yield obtained from purely hadronic rescattering without any medium effects is shown 
by dash-dotted curve. Upper right panel: The phi meson production (solid curves) and decay rates (dashed curves) for the 
process (j> <-> K~ . The thin curves are without in-medium mass modification while the thick ones include medium mass 
modification. Lower panels: cj> production (solid curves) and absorption rates (dashed curves) for different channels without 
any mass modifications; where M = ■K,p,ijj. 



curve in the upper left panel of Fig. ^, this leads to a suppressed peak in the time evolution of phi meson abundance 
but does not change much the phi meson number at later times, when compared to that in the case of including the 
initial phi mesons from string decays. 



2. rapidity distribution 

The rapidity distribution of phi mesons reconstructed from and /i^/^^ channels in central Pb-fPb collision 

at 158A GeV is shown in Fig. |^. In absence of any medium effects, the results for phi mesons from the the K^K^ 
channel (thin solid curve) are in good agreement with the NA49 experimental data (solid circles). With free-space 
decay width of phi meson, the rapidity distribution of all phi mesons from the , neglecting their scattering, is 

identical to that reconstructed from the channel (thin dashed curve). Due to rescattering or absorption of the 

kaon pairs, ~ 30% of all mesons are lost in the reconstruction, of which about 40% of the decaying kaons have at 
least one (anti)kaon that suffer elastic scattering. The maximum depletion of phi mesons from the K^K^ channel 
occurs at midrapidity where kaon-antikaon pairs undergo appreciable scattering in the dense hadronic medium. Note 
that around y = 0, phi meson from dimuon channel is about a factor 1.7 larger than from the kaon-antikaon channel, 
which is smaller than the factor of 2-4 enhancement observed in the NA50 data j|] . 
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FIG. 3. Rapidity distribution of phi meson reconstructed from K'^ K~ pairs (solid curves) and from ^~ channel (dashed 
curves) for Pb+Pb collisions at 158^ GeV at an impact parameter of & < 3.5 fm in the AMPT model. The results are for 
without (thin curves) and with (thick curves) in-medium mass modifications. The dotted curve corresponds to phi mesons 
from the dimuon channel with in-medium masses and with the phi meson number from HIJING increased by a factor of two. 
The solid circles are the NA49 experimental data |7| from the K'^ K~ channel. 



3. transverse mass spectrum 

In Fig. ^, we show the transverse mass spectra at midrapidity for phi meson from kaon-antikaon (solid squares) 
and dimuon (open squares) channel without medium effects. It is seen that at low ttit the phi meson from K~ 
channel is suppressed due to rescattering. Since the transverse momentum of a particle increases due to increasing 
number of scattering and due to pressure build-up inside the system, the decayed kaons at the early stages, which are 
predominantly scattered, thus have low transverse momenta. The transverse mass spectra can be approximately fitted 
by exp(— tot/T"). The inverse slope parameter T reported by NA50 Q from the fJ.^ channel for 1.7 < mr < 3.2 
GeV/c^ at midrapidity is T 227 ± 10 MeV, which agrees well with the AMPT model prediction of T = 228 
MeV. In contrast, the slope parameter extracted by NA49 0| from the K~^K~^ channel for 1 < mr < 2.2 GeV/c^ 
is T = 305 ±15 MeV. The depletion of reconstructed at low tjit in the AMPT model leads only to a slighter 
higher slope of T = 267 MeV, and thus is much smaller than the NA49 data. The ratio of the yields at midrapidity, 

Rimr) = Nji+K-{^cljl^4,^KK) I ^/x+/^- (r^/P^^^^), from Eqs. (||) and (0) corrected by their respective branching 
ratios, is shown in the inset of Fig. ^ The maximum suppression of ~ 40% at low rriT in the kaon channel as 
observed in the AMPT model (squares) is larger than the suppression factor of 60% found in the RQMD calculation 
This can be traced back to enhanced phi meson production from large coUisional scattering among the mesons 
leading to a peak in the (j) yield at the early stage (see Fig. ^ upper left panel). The dimuons from these phi meson 
escape the fireball freely while the K~^K^ are rescattered in the dense hadronic medium and do not contribute to the 
reconstruction of phi mesons. 



4- In-medium effect 



Modification of phi meson and kaon masses is expected to enhance the production and decay of phi meson in the 
medium and to lead to a possible further increase of the suppression factor R{mT). The abundance of phi meson as 
a function of time of the colliding system with in-medium mass modification is shown in the upper left panel of Fig. 
H (dashed curve). The abundance exhibits 20% increase at the peak, and it finally merges with the free-space value 
at large times when medium effects are negligible and chemical equilibrium in all channels sets in. This enhancement 
can be mainly attributed to the increase in the K~ —> (p production rate due to larger width r*_^^^ in the dense 
medium, as shown by thick curves in the upper right panel of Fig. 0. 
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FIG. 4. Transverse mass spectra for midrapidity (|y| < 1) phi mesons reconstructed from K~ pairs (solid symbols) and 
from jJ,'^iJ,~ channel (open symbols) for Pb+Pb collisions at 158^4 GeV at an impact parameter of 6 < 3.5 fm in the AMPT 
model. The results are for without (squares) and with (triangles) in-medium mass modifications. The solid circles are the 
NA49 data ^ for (f> — > K~ decay, and the open circles are the NA50 data ^] for (f> /i^ fi~ . In the inset is shown as 
a function of nir the ratio R{mT) for phi mesons decaying to kaon-antikaon pairs that are not scattered to those determined 
from the dimuon channel. The results are for without (squares) and with (triangles) in-medium mass modification, and with a 
further increase of phi meson number by two in the HIJING model (circles). 

The rapidity distribution of (p from kaon pairs that have not suffered any collision is shown in Fig. ^ with in- 
medium masses (thick solid curve), which is found to lie within the error bars of the NA49 data. Since free-space 
branching ratio is used to determine the phi meson abundance from the kaon channel, it is evident from Eq. that 
large in-medium width r*_^^^ leads to appreciable production and simultaneous decay of phi meson, resulting in 

dN/dy « 7.1 at y f« for all these kaon pairs. In contrast, since the branching ratio for cj) f^^fJ-^ is largely unaltered 
in the medium, the dN/dy of the reconstructed phi mesons at ?/ ~ is 4.5 and thus about 1.9 times larger than that 
from the kaon channel. 

The transverse mass spectra in the kaonic channel with medium effects (solid triangles in Fig. ^) reveals nearly 
identical slope as for the bare masses. On the other hand, the large number of dimuon production at low mx in the 
early stage of collisions with in-medium masses leads to a slightly steeper mass spectra with a slope parameter of 
T = 220 MeV. This is clearly seen in the inset of Fig. ^ (triangles) for the ratio R{m,T) where the suppression at low 
niT is ~ 33%. 

5. discussions 

Our results on the phi meson yield reconstructed from kaon and dimuon channels differs at most by about a factor 
of two in the AMPT model, which corresponds only to the lower bound to the differences found in the NA49 and 
NA50 data. For a systematic comparison with theory it is instructive to undertake an experiment where both the 
kaon and muon pairs from phi mesons are determined in the same mx and y range. If differences as large as a factor 
of four corresponding to the upper bound of the NA49 and NA50 data is indeed observed in the same experiment, this 
may then suggest that other mechanisms, such as formation of color rope or quark-gluon plasma Q, are needed for 
phi meson production. In the QGP scenario, because of copious production of ss pair, its subsequent hadronization 
may result in a dramatic increase of phi meson abundance in excess of that produced purely by hadronic rescatterings 
. The phi meson produced at the initial stage should contribute primarily to the dimuon channel. To mimic the 
effect of enhanced phi meson production in QGP, we make an ad hoc increase in the phi meson number by a factor of 
two in the HIJING model. The time evolution of A^^ is illustrated in Fig. || (dotted curve) where the phi meson yield, 
instead of exhibiting a peak, rapidly decreases to its equilibrium value which corroborates the findings of Ref. [^ . 
Including also the medium effect, the rapidity distribution from non-scattered kaon is increased by another ~ 8%, 
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whereas the yield from /i channel is considerably enhanced (dotted curve in Fig. H). The resulting suppression 
factor R{mT) at low uit is ~ 21% as is evident from the inset of Fig. ^ (circles). 

The importance of initial phi mesons to the dimuon yield is also seen in the scenario of neglecting phi mesons from 
the initial string decays. In this case, we find that the phi number obtained from the K^K^ channel is essentially 
unaffected due to chemical equilibration via the hadronic scatterings, but the phi meson reconstructed from dimuons 
is, however, significantly suppressed. 

B. RHIC 
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FIG. 5. The rapidity distribution (top panel) and the transverse mass spectra (bottom panel) for midrapidity {\y\ < 0.5) 
phi mesons reconstructed from K'^ K~ pairs (solid curves) and from iJ.'^fi~ channel (dashed curves) for Au+Au collisions at 
RHIC energy of -^s — 130A GeV at an impact parameter of 6 < 5.3 fm in the AMPT model. The solid circles are the STAR 
experimental data for — 11% central collisions for cj> reconstructed from K~ decay. 

In the AMPT model, heavy ion collisions at RHIC differ from that at SPS mainly in the increasing effect of minijet 
gluons. However, the initial heavy ion collision dynamics at RHIC remains dominated by soft string fragmentation 
processes as at SPS. In Fig. || (top panel), we present results from the AMPT model for the rapidity distribution of phi 
mesons reconstructed from K'^K~ and jJL^ ^~ pairs for central Au+Au collisions at the RHIC energy of -y/s = 130A 
GeV. Of the total phi meson yield of 37 per event, about 24% of these are lost in the kaonic channel by hadronic 
rescattering and absorption. Note that in spite of enhanced phi meson production in the early stages of the collision 
as compared to the SPS energy, its abundance at ?/ ~ is only about a factor of 1.5 large in the dimuon channel. 
This may be attributed to fewer kaons lost by rescattering at RHIC as they can escape rapidly out of the collision 
zone unperturbed due to of their large energies. It is seen that phi meson multiplicities reconstructed from K^K^ 
pairs at midrapidity in the AMPT model is consistent with the STAR data p| . 

The transverse mass spectra of 4> meson from the two channels are shown in the bottom panel of Fig. |[ and 
compared with the STAR data |^l[] for the K^K~ channel. Compared to a slope parameter of T = 379 ± 50 
MeV in the data, the AMPT model predicts a smaller value of T = 335 MeV in the kaonic channel in the range 
< Tot — < 1 GeV. The slope parameter for phi mesons determined from the /i^/i^ channel is T = 297 MeV 
resulting in a suppression factor at low tot of R{mT) — 58%. Note that at both SPS and RHIC energy, the slope 
parameter of phi mesons from K~ decay is smaller compared to the experimental data. We have increased the 
phi meson elastic cross sections with baryons and mesons to 8.3 mb, which corresponds to the upper bound estimated 
from phi meson photoproduction data ||2l| , and find that this increases only slightly the slope parameter. On the 
other hand, an enhanced flow is expected to be generated by converting the strings produced from soft processes 
into interacting partons in the high energy density created in ultra-relativistic heavy ion collisions [p^ . Since in the 
default HIJING model the strings are basically noninteracting, the latter modification was found to reproduce the 
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large elliptic flow observed in Au+Au collisions at -Js = 130 A GeV [|T^. Results based on this new AMPT model 
will be reported in the future. 

V. SUMMARY AND CONCLUSIONS 

In summary, we have investigated in a multiphase transport (AMPT) model the production of phi mesons recon- 
structed from K'^ K~ and fJ-'^ n~ decays. Considering all possible collisions of cf) with the hadronic medium, we find 
the phi meson yield from dimuon channel is about a factor of 1.7 higher than the kaon-antikaon channel for central 
Pb+Pb collisions at the SPS energy. The kaons originating from phi meson decay in the dense hadronic medium are 
predominantly scattered and absorbed and thus do not contribute to the phi meson yield, whereas all the ^i" pair 
freely escape. Inclusion of in-medium modification of masses of kaons and phi mesons enhances the phi meson width 
and results in a factor of two increase of the phi meson yield from the dimuon channel compared to the kaon-antikaon 
channel. However, strikingly different results were found by two experiments: The NA50 [^j measured cj) /^^M^ 
and found its yield a factor of 2-4 larger with a smaller slope parameter compared to ^ K^K^ measurements 
by NA49 0] for central Pb-|-Pb collisions at 158A GeV. Such a large enhancement of four at low transverse mass 
from dimuon over the kaon channel could only be explained in the AMPT model by an ad hoc increase of phi meson 
number by about three from the initial stage. This may be suggestive of other mechanisms, such as the formation 
of color ropes or quark- gluon plasma in these collisions. To pin down these differences an experiment measuring phi 
meson from both dimuon and kaon decays in the same kinematic range is therefore called for. For heavy ion collisions 
at RHIC energy of y/s = 13QA GeV, the AMPT model gives the yield and slope parameter of phi mesons from the 
K^K^ channel that are in reasonable agreement with the STAR data. Similar to the results for heavy ion collisions 
at SPS, the phi meson yield from the dimuon channel is about 1.5 times larger than that from the K^K~ channel 
in the absence of additional phi meson production due to the formation of the initial partonic stage. This number 
is expected to increase significantly if we include such enhanced production of initial phi mesons. Indeed at RHIC, 
large energy densities are expected to be reached, and a partonic matter with a longer lifetime and occupying a larger 
volume could be formed. Compared to SPS energy, larger ss production in the partonic stage and its subsequent 
hadronization would therefore result in an enhanced phi meson abundance, especially at the early stages. These phi 
mesons should contribute dominantly in the dimuon channel while most of them will escape detection in the kaonic 
channel. It is thus of great interest to have experimental data on phi mesons from the dimuon measurement at the 
RHIC energy. 
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